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PREFACE 


The objectives of this study were to determine the feasibility of 
using data from the Electrically Scanning Microwave Radiometer (ESMR) to 
1) characterize the iroisture condition of the upper soil zone into one 
of two or three broad classes* and 2) delineate regions of snow cover. 

With respect to the first of these, comparisons of ESWl-5 bright- 
ness temperatures to soil moisture estinates conputed by a daily bud- 
geting scheme »«re made. These correspond to locations in the U.S. 

Great Plains in the August-October time period. Despite the drawbacks 
associated with the time period (relatively high vegetative cover and 
few substantial rain events), feasibility under a limited set of conditions 
was denonstrated. The major limitation is the vegetative cover, which 
must be minimal. 

Study of snow pack delineation was limited to a single case, and 
deironstrated feasibility for ESMR-6. 


iii 


TABLE OF CONTENTS 


ABSTRACT 

PREFACE 

INTRODUCTION 

BACKGROUND AND APPROACH 

SOIL MOISTURE BUDGET 

PROCESSING OF ESMR DATA 

ANALYSIS OF SOIL MOISTURE 

SUMMARY OF SOIL MOISTURE ANALYSIS 

SNOW COVER DELINEATION 

CONCLUSIONS 

RECOMMENDATION 


REFERENCES 



INTRODUCTION 


Recent investigations have demonstrated that data from a satellite- 
borne radiometer, sensitive to passive radiation at microwave frequencies, 
can provide Information of great potential value to such diverse Inter- 
ests as ship operations, agriculture, and hydrology. In the earlier 
Investigations, the radlaneters were flowi on aircraft. The first 
satellite-borne microwave radiometer was the Electrically Scanning 
Microwave Radiometer (ESMR) carried onboard the Nimbus-5 satellite. 

ESMl-5 is sensitive to radiation In a narrow band centered at 19.35 ^2 
(wavelength 1.55 cm). At this frequency, the variations In the observed 
radiation are strongly Influenced by variations in the emlssivity of the 
surface. The major contribution to emissivity is liquid water with its 
extren«ly large dielectric constant at microwave frequencies. Many 
investigations have been concerned with surface phenottwna over the 
ocean. Sea-surface wind speed can be inferred at microwave frequencies 
through the raising of emitted energy due mainly to wind generated foam. 
Investigations of this phenorrrenon have been conducted by Williams 
(1969), Nordberg et al. (1971), and Sabatini (1975), to naire a few. 
Boundaries of large-scale fields of sea Ice can be mapped, even In the 
presence of clouds, as has been shown by Gloersen et al. (1973) and 
Sabatini et al. (1975). Sabatini et al. (1975) have also shown that 
ESMR-5 data is useful In delineating areas of rainfall over the ocean. 
Detection of rainfall over land Is desirable from agricultural and 
hydrological considerations, but Is more difficult due to the highly 
variable nature f the land background. Wilheit et al. (1975) show that 
under certain conaltlons rainfall rate can be roughly estimated over a 
land surface. Meneely (1975) showed that in relatively level regions of 
minimal vegetative cover, areas of recent rainfall could be delineated 
and relative rainfall amounts estimated using ESMR-5 data. Such de- 
li neaHorTFesu Its from the lowering of the surface emissivity as water 
permeates the upper portion of the soil. That rainfall study comple- 
mented investigations of soil moisture detection with airborne radio- 
meters by Schmugge et al. (1974 and 1976). 

A knowledge of the moisture content of the upper soil zone is 
potentially useful in another context, namely as a key in assessing 
large-scale moisture flux into the atmosphere, which could serve as an 
Input to a n^teorologlcal atnwspheric model. It was this potential for 
determining large-scale moisture flux which provided the motivation for 
the study reported here. 

As specified in the RFP for this study, data from the Nimbus-6 
ESMR, due for launch In early 1975, would be used to examine two aspects 
of ground surface condition. The major emphasis would be on soil mois- 
ture, with a brief look at snow cover detection. The test period for 
the soil moisture study was to be 1 May to 30 June 1975 and for snow 
cover 15 December 1975 to February 15 1976. 

The key element in the study was the calculation of soil moisture 
values, to be used In a "ground truth" context. In situ measurements of 
soil moisture are difficult to take and to calibrate. A soil moisture 
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budgeting schane previously developed and thoroughly tested by Earth 
Satellite Corporation m& to be used to compute the water content of the 
upper portion of the soil on a daily basis. The area to be studied was 
in the Great Plains, extending from southern Canada to northern Mexico. 

The Nimbus-6 ESMR differs from ESMR-5 in three important aspects: 
a) it is sensitive at a shorter wavelength, 0.8 cmi b) it scans along a 
conical arc so it makes an essentially constant scan angle with the 
surface! and c) it senses the horizontally and vertically polarized 
components of the radiation. The actual launch of Nimbus-6 was delayed 
for several months, precluding the originally planned soil moisture 
study period. Early problems with availability of properly calibrated 
ESMR-6 data led to a decision by NASA to use data from ESWl-5. Accord- 
ingly, ESMR-5 data for the periods September through October 1974 and 
August through October 1975 were supplied by NASA for the soil moisture 
analysis. The resulting analysis demonstrated the limitations imposed 
by the presence of a crop canopy and lack of many significant rainfalls 
at this time of year. 
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BACKGROUND AND APPROACH 


The Intensity of radiation emitted by a surface Is proportional to 
the product of the emissivlty of the surface and its temperature. This 
product Is coiwiwnly called brightness temperature (To). For satellite- 
borne radiometers, the received radiation Is modified somewhat by atten- 
uation and emission within the atmosphere. For non-precipitating atmos- 
pheres, this effect Is still quite small and surface emissivlty remains 
the major contribution to variations in the sensed Tg, in the type of 
application rm>t1vat1ng the current study, narrely specification of Inputs 
to an atmospheric model, an average value for Tg would be determined for 
each of a nunter of large computational cells. As originally conceived 
for this study, each cell would encompass 2.5“ latitude by 2.5® longi- 
tude (about 50,000 square klloneters). It would therefore contain soite 
80 ESMR-5 spots (assuming the sub-satellite spot size of 25 x 25 km). 
Fifty-six of these computational cells covered the area of interest for 
this study. 

The daily soil moisture was to be calculated based on a scheme 
originally developed by Baler and Robertson (1966) 1n which the soil 
profile is divided into several zones. Moisture is added to or subtracted 
from each zone based on precipitation, atmospheric demand, soil character- 
istics, estimated crop root development, and the amount of water already 
in the zone. For the study reported here, it Is the uppermost zone that 
Is of fundamental interest. The scheme is discussed more fully beginning 
on page 5. One of the basic inputs for this computation 1s the dally 
atmospheric demand, or potential evapotranspi ration (ETP). Effective 
definition of this parameter requires rreteorological data at six-hour 
intervals or less. The density of stations reporting data in the required 
format was such that most of the computational cells would have had only 
one or two stations within them. It was beyond the scope of this study 
to perform the operations necessary to interpolate the data to the fine 
grid mesh needed to accurately characterize each cell. This is particu- 
larly true with respect to rainfall which is highly spatially variable. 
Simple sampling theory such as Tchebycheff 's inequality readily shows 
that such a cell cannot be accurately characterized using only one or 
two points. The following numerical example Illustrating this fact was 
taken from some calculations performed for the northern Great Plains 
during 1975. In these calculations, the weather paran^ters had been 
determined for grid points with spacing of 12.5 miles of which there are 
about 110 to 120 in a 2.5° by 2.5° cell. Data for nine such 2.5“ cells 
in the region were examined for the period August 24-30. For the 63 
cases (nine cells, seven days), the average moisture content of the 
uppermost 5 percent of the soil profile varied from as low as 0.6mm to 
as high as 8.2mn, or nearly equal to the assumed saturation value of 
8.75mm. The standard deviation within the cells ranged from 0.6 to 
3.1mm. Application of Tchebycheff 's inequality showed that specification of 
the mean cell top zone moisture to within 1.0mm will be achieved at 
least 90 percent of the time with a sample of that many (110-120) points. 

For a standard deviation of 2.0, accuracy to within l.Oran 50 percent of the 
time still requires eight points within the cell. 
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In addition to these spatial variations due to rainfall patterns, 
there are the effects of soil characteristics, type of terrain, land 
use, and amount of vegetative cover. All these factors Interact to 
produce great variations In the Tg-molsture relationship, making It 
clearly Inappropriate to conduct this soil moisture Investigation based 
on the 2.5* cells. The moisture budget calculations were therefore 
performed for points corresponding to the actual locations of all weather 
stations reporting the required n^teorologlcal data. The resulting 
moisture values were compared to spatially averaged Tg values from ESMR- 
5 In the Immediate vicinity of the station. (At least four values 
surrounding each station location were used to minimize the effect of 
Instrument noise.) As reported by Schmugge et al. (1976), research has 
shown that microwave detection of soil n»1sture at shorter wavelengths 
Is strongly affected by the presence of a vegetative cover. This was 
clearly Indicated In the detection of rainfall over land areas by ESMR-5 
as reported by Meneely (1975). That study examined two heavy rain 
events which occurred In the most Intensively agricultural portions of 
Illinois and Indiana In dune and duly of 1973. Detailed soil moisture 
calculations using the Baler and Robertson model (see pages 5-7) were 
performed at six locations. Data from nine ESMR-5 passes over the area 
were used. A growth model for corn was used In the analysis to estimate 
the vegetative cover. A summary of the results appears In Figure 1 and 
shows the observed brightness temperatures plotted against the frac- 
tional saturation of the uppermost soil zone. The open symbols corres- 
pond to the early dune event, for which there was little vegetative 
cover (a fact confirmed by LANDSAT-I Imagery). A strong Tg decrease 
with Increasing moisture was noted, averaging about 40 K from dry to 
saturated (field capacity). The filled symbols represent the duly case, 
for which the vegetative cover was near 100 percent. Here the Tg 
decrease was only about 17 K. 

It was decided for the study reported here that imagery from 
LANDSAT-I and II would be used where readily available to define vege- 
tative cover within 25 km of each station. Figure 2 is a map showing 
the location of the 52 stations for which soil moisture calculations 
were made. The 42 stations for which LANDSAT images were examined are 
indicated by circles. Figure 3 presents the estimated percent of vege- 
tative cover at each of these stations. In the southern portions of the 
study region, imagery for the desired August to October time period was 
not so readily available, so late springtime imagery was substituted. 
These values are enclosed in parentheses. The analysis on Figure 3 
shows a general geographic consistency and an expected east to west 
decrease in the cover. 
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SOIL MOISTURE BUDGET 


The scheim used to calculate the dally soil tm)1sture profile was 
developed by Baler and Robertson (1966) and designated "Versatile Bud- 
get” (VB) by them. This sch^ divided the total moisture holding 
capacity of the soil Into several zones relating to a wheat plant root 
structure. Water was subtracted from different depths each day through 
coefficients which related plant rooting characteristics and soil v^ter 
evaporation to the dally atmospheric demand. It was the top-most zone 
that was of Interest to this study. This zone was assumed to contain 
five percent of the total water holding capacity (wilt point to field 
capacity) of the soil profile. For a typical wheat growing soil of one 
meter depth, this represented a physical zone depth of about five centi- 
meters. 

The VB model for calculating the actual dally nwlsture extraction 
from a given soil zone was given by 


ET, • K.(^).Z.ETP, 



-1 



where 



S 

Actual evapotranspiratlon for the layer on day 1 
ending at sunrise on day 1+1. 


K 

s 

Crop/soil coefficient. 


Sr' 

3 

Available soil moisture in the zone at the beginning 
of day 1-1. 


c 

= 

Capacity for pi ant- avail able water In the zone. 


z 

= 

Adjustment factor for soil drying characteristics. 


FTP. 

s 

Evapotranspiratlon potential (atmospheric demand) 
for day 1. 


The net soil nwlsture In a zone at the end of day 1 was given by 


s, • s,., - ET, + I 


where 

I = The Infiltration of rainfall Into the zone. If the 
total rainfall for the day was less than 25.4nffi, 
then the total was assun^d to be available for 
infiltration. Otherwise, the amount available for 
infiltration was given by 




fjL) 

r’-') 

125.4' 

\ %* f 


R, » Rainfall in mm for 24 hours ending the morning of 
$• ’ day 1+1. 

1-1 s Fraction of soil moisture capacity available in top 
S* zone at end of day 1-1. 

1a was used to bring the top zone up to capacity. If any excess 
remained, it was used to fill the lower zones In turn. Any remaining 
amount was added to run-off. 


The crop/soil coefficient (K.) was taken to be 0.40 for the top 
zone. This was representative of most crops in the late growing season 
or of bare ground, and implied that 40 percent of the atnwspheric demand 
(ETP) would act on the top zone. K represents the average coefficient 
for the zone. At the top surface K is nearly 1.0, but decreases rapidly 
with depth. The soil crop-available moisture capacity (C) was a function 
of the soil type and referred to the amount of nwisture in the range 
fr«n field capacity (-1/3 bar) to permanent wilting (-15 bar). Two 
general types of soils were considered for this study. One was a clay 
loam or silt loam with a top zone capacity of 8.75mm. The other was a 
sandy soil with a top zone capacity of 5.75mm. The adjustment factor 
for soil drying characteristics (Z), commonly called the "dry-down 
curve" described the ability of the soil te give up its plant-available 
water. There are numerous viewpoints on how Z varies with the amount of 
plant-available water. The relationship used must generally be based on 
past experience. The following relationships were used In this study. 

In sandy soils, all the water was assumed to be freely accessible to 
meet the atnraspheric demand. For clay or silt loams, three dry-down 
curves were defined and computations made for each. Curve "G" assuTOd 
that the fraction of the water content available to n«et the atmospheric 
demand decreased linearly from 1.0 when the zone was saturated to zero 
when the zone was totally dry. Curve "E" assuired that all water was 
available until 35 percent of zone capacity remained. Then the avail- 
ability ratio decreased sharply and logarithmically to zero at zero 
content. Curve "E-normal" was similar to "E" except the transition 
ratio was 25 percent instead of 35 percent. 

The evapotranspiration potential (ETP) was calculated using the 
n^thod of Penman (1948) 


ETP ^ 


A + <T 
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where 


A 

hu 

R 




et 


w 


Slope of the saturation vapor pressure curve (mb/'^K). 
Heat of vaporization of water. „ 

Net radiation for the day (cal/cm^). 

Psychrometric constant. 

Wind effect on evaporation. 


In order to more accurately assess the effects of temporal var- 
iations in meteorological parameters, the day was divided into four 
periods of six hours, each centered about a meteorological observation. 
The parameters used were air tempe'^ature, dew point temperature, wind 
speed, and cloud types and amounts. 


The net radiation consisted of two components: net incoming solar 

and net outgoing longwave. The solar component was determined by inte- 
grating using half-hour time steps to account for the earth's rotation. 
Attenuation and scattering by dust, water vapor, and clouds, and reflec- 
tion from the surface were included. The longwave component was sub- 
tractive and was a function of the surface temperature (assumed to equal 
the air temperature) and the water vapor content of the air. Reflection 
by clouds was also included. The wind effect was a function of the wind 
speed and the difference between the actual and saturation vapor pres- 
sures. 

The meteorological data used in these calculations were obtained on 
magnetic tape from the National Climatic Center in Asheville, North 
Carol ina. 
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PROCtSSING OF ESMR DATA 






Uncertainties in satellite attitude make wholly automated qeo- 
referencing of Tg spot values unacceptable. A technique has evolved in 
which tnis is done manually using q~idded overlays developed T. 'm the 
orbital geometry, ’’’he ESMR calibrated brightness temperature tapes are 
read and the Tg spot values ure printed in a format which produces a 
proportionally correct relationship between the spacing of the spots 
within a scan line and the spacing of the scan lines themselves. De- 
pending on the sensor (ESMR-5 or 6), the hemisphere (northern or southern), 
and the direction of satellite travel (northbound or southbound), the 
proper overlay is selected and initially positioned using the printed 
ephemeris data for each scan. Water bodies and coastlines, which are 
readily visible in the Tg field, are outlined and the overlay is shifted 
as needed to place these at the proper coordinates. The latitude- 
longitude grid is then transcribed onto the Tg map with an error gener- 
ally much less than the spot spacing. 

Each ESMR-5 scan line consists of 78 individual spots. The quality 
of the data is seriously degraded at the ends of the scan due to in- 
creased atmospheric contribution and poorer resolution. For this 
reason, only spots which were more than ten positions from either end of 
the scan were used. The available ESMR-5 orbits were screened and those 
which included at least a part of the study area within the acceptable 
region were selected for further analysis. A total of 21 daytime and 11 
nighttime passes was selected. 
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ANALYSIS OF SOIL MOISTURE 


PAGU liLAMv NuT FILM .:iD 


The analysis of daytime and nighttime data could not be combined 
because of the large differences in surface temperature and the diffi- 
culty in unambiguously removing their effects. Only the analysis of 
daytime data is reported here; the nighttime data set was too small to 
permit significant conclusions. The soil moisture budget calculations 
at each station were started a number of days prior to the first ESMR-5 
pass (19 days prior in 1974, 25 days in 1975). This allowed the effect 
of uncertainties in the initial estimates of soil moisture to diminish. 
Determination of the proper soil type and dry-down curve for each station 
is a difficult task at best. For this reason, all four combinations of 
soil type and dry-down curve were used at each station. It was hoped 
that the resulting comparisons would indicate which soil description was 
most appropriate at each station. It was found, however, that the 
results were relatively insensitive to the soil parameters except on the 
one or two days following rainfall. There were too few instances of 
significant rainfall with ESMR-5 data available on immediately following 
days to permit this assessment. Table 1 is an example of the output of 
the computation for Williston, North Dakota. The soil moisture in the 
top zone was expressed as a fraction of the capacity. To compress the 
output to one page, only those days on which or near which a usable 
ESMR-5 pass occurred were printed. The hand-tabulated brightness tem- 
perature values represented the average of the four spot values surrounding 
each station. Orbits not tabulated indicate that the station was not 
within an acceptable portion of the scan. The only available daily 
rainfall totals were for the 24-hour period ending at midnight. Since 
the daytime ESMR-5 passes occurred near noon, it was not appropriate to 
use data from a day on which rain fell between midnight and noon. 

Station reports of rainfall during those hours were indicated by x's. 

Data from those days were not used (i.e., October 5, 1974 at Williston). 

For each station and for each soil/dry-down combination, simple 
plots of Tg versus fractional saturation in the top zone (f^) were 
generated and correlation coefficients were calculated. Figures 4 
through 8 present examples for five of the stations. The vertical axes 
are Tg and the horizontal are the ratio between the actual moisture 
content of the top soil zone and its field capacity. The number printed 
at each plotted location is the number of data points occurring at that 
location. It was felt that much of the scatter appearing in the oata 
was attributable to day-by-day changes in the soil temperature. No soil 
temperature data were available; however, an attempt was made to com- 
pensate for its effects by assuming that the mid-day soil surface temp- 
erature could be approximated by the daily maximum air temperature. 

These results are presented in the lower tier of plots for each station. 
While this did reduce the scatter for some stations, it increased it for 
others, and the net result for all stations aggregated was to decrease 
the magnitude of the correlation coefficient somev/hat. A more .'-romisinq 
approach might estimate soil temperature based on calculations of the 
solar energy reaching the soil up to the time of satellite passage. 

Kansas City, Missouri (Figure 4) had the highest correlation of all the 
stations, particularly when the temperature adjustment was included. 

The highest correlation was obtained for dry-down curve 6, the lowest 
for sandy soil. Correlation coefficients were also produced by ex- 
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eluding the driest data points. This was done because at some local- 
ities there seemed to be too many points grouped at this end of the 
curve to allow a n^aningful result. In nrost cases, however, the corre- 
lation did not change much and the reduced sample contained too few data 
points to be rr^aningful. Des Moines, Iowa (Figure 5) is another station 
which displayed a relatively high correlation coefficient, but in this 
case the lowest correlation was for dry-down curve G. While these two 
stations had fairly high correlation coefficients, the sensitivity of Tg 
to iroisture content was not great - about 15 K from dry to saturated. 

This Is likely due to the relatively high vegetative cover (estimated at 
50 percent and 65 percent, respectively). Figure 6 presents another 
station with a high correlation coefficient - Williston, North Dakota. 

In this case, there is a much greater sensitivity to moisture content - 
about 30 K from dry to saturated. This is consistent with the estimated 
vegetative cover of about 25 percent. Two other stations, only about 
100 miles apart, are shown to denwnstratc the gross differences in 
response which they display. One of the stations is Oklahoma City, 
Oklahoma, which displays a high correlation coefficient and a dry to wet 
range of about 25 K; the other is Tulsa, Oklahoma, which displays no 
correlation at all. The fraction of vegetative cover is similar at 
these two locations, but the discrepancy probably arises fran the fact 
that the Tulsa area is highly forested, while the Oklahoma City area Is 
mostly range and cropland. 

Figure 9 is a map of the test region on which the correlation 
'^efficients have been plotted and analyzed. Several stations across 
the center of the region had very little rainfall during the test per- 
iod. In these cases, insufficient data was available at higher moisture 
levels to allow rr^aningful results. These stations are indicated by 
asterisks. There seems to be a tendency for the regions where grain 
crops are predominant to be the most highly correlated. This includes 
the spring wheat belt from north central and eastern Montana across 
North Dakota to western Minnesota, the western portions of the corn 
belt, and the winter wheat belt through southern Kansas, western Oklahoma, 
and the Texas panhandle. In the time period being considered, the 
spring wheat will have been mostly harvested and the winter wheat will 
be in an early stage of growth, implying a fair percentage of bare soil. 

In the corn belt, harvesting will be at an early stage; however, in the 
western portions there may well be a great enough mix of crops to 
assure some bare soil. 


Figure 10 presents a plot of all the Tn versus fractional sat- 
uration data points for the 52 stations. There is a great amount of 
scatter in the data and only a weakly discernible slope. Figures 11 
through 13 present similar plots for the three regions of highest corre- 
lation coefficients visible in Figure 9. In these cases, the trends in 
Tg are much more clearly discernible. Many of the outlying points, 
particularly those with high Tg values at higher moisture levels are not 
present in Figures 11 through 13. The dashed line on Figure 10 (all 
points) defines the limit above which no points exist on the three 
following plots. These outlying points represent, for the most part, 
locations with heavier natural vegetation lying outside the major agri- 
cultural regions. These grasslands or forests mask the return from the 
wet soil benciith. In the spring wheat and winter wheat regions (Figures 11 
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and 13), the To changes fr«n wilt point to field capacity are 24.5 K and 
23.3 K, respectively as determined by least-squares regressions. These 
correspond to average vegetative cover estimates of about 40 percent. 

In the corn belt region {Figure 12) where the average vegetative cover 
is estimated at 55 percent, the Tg decreases by only 11.7 K. 

In these three figures, the most notable outlying points are those 
associated with low calculated soil moisture in conjunction with low 
observed brightness temperature. These are likely manifestations of 
shortcomings in the soil moisture model. In particular, the eight 
points lying belw# the dashed line of Figure 11 <p«re subjected to fur- 
ther scrutiny. 

The follwing is a list of the locations, dates, observed noon-tiire 
Tb» and fractional saturation (fg) estin»tes for the upper soil zone: 


Station 

Date (1975) 



Minot, NO 

Oct 1 

247 

0.22 

Williston, NO 

Oct 5 

252 

0.34 

Billings. MT 

Oct 19 

253 

0.13 

Havre, MT 

Oct 19 

252 

0.17 

Glasgow, MT 

Oct 21 

245 

0.19 

Havre, MT 

Oct 21 

251 

0.12 

Minot, NO 

Oct 21 

243 

0.22 

Williston, NO 

Oct 21 

243 

0.41 


These observations have two things in common. First, they all occur 
in October 1975 and, second, all of the stations are located in the 
northern Great Plains. This proximity in space and time suggests that 
the anomalies are a result of real physical causes and not merely random 
deviations from the average. 

Addressing the Williston, Oct. 5 observation first, an inspection 
of the rainfall data indicates that a light rainfall (2mm) occurred on 
the evening of Oct. 4. The soil moisture model estimated that this 
rainfall brought the level of saturation to 0.34 in the top zone. The 
design of the model is such that the moisture is introduced uniformly 
throughout the (approximately) five centimeter depth of the zone. 

However, since the rainfall occurred just the evening before the satellite 
observation, the moisture could reasonably be expected to be concentrated 
in the top few millimeters of the soil. In such a case, the satellite 
would detect the wet soil surface and would produce readings of Tg which 
correspond to greater soil moisture than the model estimates. 

The above explanation might be applied to the Minot, Oct 1 obser- 
vation. A trace of rain was detected at Minot on Sept. 30. This could 
indicate that the Minot area had a light rain which failed to infiltrate 
or evaporate significantly by the noon observation on October 1. The 
morning of that day was mostly cloudy and cool with high relative 
humidity, further supporting this assumption. 


An inspection of the meteorological data for the area for the week 
prior to Oct 21 indicates a cause different than the one above for the 
remaining six anomalous observation points. For a six day period, 

Oct. 11-16, precipitation occurred in North Dakota and Montana. The 
total amounts for the stations in the list varied from slightly over 
12mm at Billings to over 40mm at Minot. Very importantly, a sizeable 
portion of this precipitation was snowfall. Although warm temperature 
caused rapid melting and none of the stations reported any snow cover 
after Oct. 17, the snowfall could cause the observed Tg to be lower than 
otherwise expected in two ways. First, patches of snow might still 
exist at observation time in the areas near the stations. The patches 
of snow would contaminate the Tg observation because of the cold tem- 
perature and different emissivity of the snow surface. Second, and more 
importantly, the snow cover acts as a source of moisture until it is 
melted. The soil moisture model makes no allowance for this effect of 
snow cover and begins depletion immediately after the last precipi- 
tation, using the energy which melts the snow to evapo>"ate the water 
instead. This premature depletion of soil moisture by the model when 
precipitation falls as snow will cause under-estimation of the amount of 
soil moisture for several days after the last precipitation and can 
result in such anomalous points as were observed on Oct. 19 and Oct. 21. 


15 



SUMMARY OF SOIL MOISTURE ANALYSIS 
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At short wavelengths, such as the 1.55 cr value for ESMR-5, the 
absence of a vegetation canopy is critical to its ability to detect soil 
moisture variations. This is clearly illustrated by Schmugge et al. 
(1976) in comparing ESMR-5 Tg values with those for two radiometers 
carried on SLylab. The Skylab radiometers had wavelengths of 2.2 cm and 
21.0 cm. The data were for a 300 km swath across north-central Texas on 
June 5, 1973. The soil moisture in the top 2.5 cm varied from 15 percent 
to 70 percent of field capacity. The 21 cm radiometer displayed a Tg 
range of 45 K, while the range at 2.2 cm was 15 K. The ESMR-5 range 
was no more than 5 K. 

The cases presented here showed that only in the presence of mini- 
mal vegetative cover could ESMR-5 meet the basic objective motivating 
this study. That was the ability to specify whether the soil moisture 
was at or below the wilt point (f^ <_ 0), saturated (f^ ^1) or between 
these extremes. If the slope of the Tg variation is small with respect 
to the scatter in the data, then little confidence will exist in the 
assignment of a moisture class based on an observed Tg. Even in those 
cases of lesser vegetative cover (Figures 1-open symbols, 11 and 13) 
there still exists a great deal of scatter. There are a great many 
possible sources for this scatter, including: 

a) variations in ground temperature among the data points; 

b) variations in the water content of the atmosphere; 

c) uncertainties associated with the soil moisture model, in 
particular, the large thickness (up to 5 cm) of the uppermost 
zone in comparison to the wavelength of ESMR-5 (1.55 cm) and 
the characterization of the atmospheric demand response ot 
the entire zone by a single coefficient; 

d) local topographical effects (as they influence water retention 
and runoff). 
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SNOW COVER DELINEATION 

Because the microwave emissivity of a snow cover can differ mar- 
kedly under many conditions from nearby bare ground, a microwave radio- 
meter provides the opportunity to map the extent of a snow pack, even in 
the presence of clouds. A brief analysis of some ESMR-6 data was con- 
ducted to verify that useful information can indeed be obtained. The 
orbit selected was a daytime pass over the central U.S. on December 9, 
1975. Figure 14 is a portion of the analyzed map of the horizontally- 
polarized Tg (Thor) centered over east-central Minnesota. The outline 
of a portion of Lake Superior was indicated by the 200 K isotherm. 
Circled numbers represent snow depths in centimeters as reported at 
cooperative observing stations in the NOAA publication Cl imatological 
Data . This snow pack was several days old and some prior meeting had 
occurred, although air temperatures were well oelow freezing on the day 
of this image. A sharp change in T^qr was noted at the southern boun- 
dary of the snow field, with snow-covered areas being 15 K to 30 K 
lower. The 230 K isotherm (shaded) appears to be a good delineator of 
this boundary, where the snow water equivalent is highest and where the 
greatest melting has occurred. The 230 K line to the north does not 
indicate a boundary, but merely reflects a northward increase in the 
pack emissivity as its surface characteristics change. 

Analysis of the vertically polarized component showed similar 
trends but the range of magnitudes was less. Analysis of the difference 
of the two Td components demonstrated little relationship to the snow 
field. 
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CONCLUSION S 

8.1 Soil Moisture 


The most important conclusion resulting from this study is: 

A microwave radiometer of short wavelength (less ^han 
2.0 cm), such as ESMR-5, appears to be able to provide useful 
information on the moisture condition of the upper soil zone 
provided that all of the following conditions are met: 

1. Little or no vegetative cover (less than ^0 percent) 
must exist over a wide region. 

This implies that the areas must be either 
arid or heavily agricultural with either crop- 
land or mixed crop and range land, and confined 
to the time of year whe?i fields are bare or 
nearly bare. 

2. Soil surface temperature must be estimated within 
the region. 

Since the soil moisture directly affects 
the emissivity, those influences on Tg which 
are not manifested in emissivity must be accounted 
for. Soil temperature is a major item in this 
category. 

3. A regional Tg/moisture relationship calibration must 
be made. 


This will allow for some of the effects of 
differing soil characteristics, drainage effi- 
ciency, and land use mix to be factored out. 

Much more successful large-area soil moisture determinations 
over a greater part of the year v/ould be possible with a radiometer 
sensitive at a considerably greater wavelength. 

Determination of the feasibility of extrapolation into the 
winter season or into arid regions with wide mixes of terrain type 
requires further study. 


8 . 2 Snow Mapping 

Mapping the active (southern; margin of a snow pack using 
microwave radiometers is feasible. A horizontally-polarized Tg 
component of 230 K seems to be a reasonable delimiter. 
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RECOMMENDATION 

The overriding recommendation resulting from this study is for a 
more detailed evaluation based on data for the months March through June 
in the northern Great Plains and the Midwest. In such a study, morning 
solar radiation calculations would allow reasonable estimates of soil 
surface temperature to be made. Three or four regions would be defined 
and the data aggregated to that level to permit assessment of regional 
differences. 
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BRIGHTNESS TEMPERATURE 



0 .1 .2 .3 .4 .6 .6 .7 .8 .9 1.0 

FRACTIONAL SATURATION OF TOP ZONE 
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FIGURE 1 : Brightness temperature versus fractional saturation of top zone for 

six locations in Illinois and Indiana, June ■ July 1973. Open symbols 
indicate less than 20% crop canopy. Filled symbols indicate more 
than 70% crop canopy. 
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FIGURE 2: Locations of the 52 stations for which soil moisture calculations were made are 
indicated by X. The 42 stations for which LANDSAT imagery was obtained 
are circled. 
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FIGURE 3; Estimated percentage of vegetative cbver based on LAIMDSAT imagery. 
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FIGURE 4 Da-t.me b..gfMne« temf^rature 'Tr. « K) vs. fract-onal saturation in top soil zone for Kansas City, Missouri. Estimated vegetative cover is 50^.. 
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FIGURE 5 Daytime brightness temperature (fp, ® K) vs. fractional saturation tn top soil rone for Des Moines. Iowa. Estimated vegetative cover is 65%. 
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FIGURE 6; Daytime brightness temperatun {Tg, “ K) vs. fractional saturation in top soil zone lor Williston. N. Dakota. Estimated vegetative cover is 25%. 
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FIGURE 7: Uayt.me brightness temperature <Tb, « K) vs. fractional saturation m top soil zone for Oklahoma City, Oklahoma. Estimated vegetative cover is 55% 
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FIGURE 8: Daytime brightness temperature (Tg, ° K) vs. fractional saturation in top soil zone for Tulsa, Oklahoma. Estimated vegetative cover is 60%. 
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FIGURE 9: Map of correlation coefficients of brightness temperature vs. fractional saturation 
of soil top zone. Asterisks f ♦) indicate those stations with too little rainfall to 
permit computation. 
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FIGURE 10. Daytime brightness temperature versus fractional saturation of top zone 
for 52 stations in the United States Great Plains for September - October 
1974 and August - October 1975. Numerals give number of data points 
at each position. Asterisk (*) indicates ten or more.Dashed line defines 

upper limit of data in Figures 11 through 13. 
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FRACTIONAL SATURATION OF TOP ZONE 


FIGURE 11. Daytime brightness temperature versus fractional saturation of top zone 
for eight stations in the United States spring wheat belt for September - 
October 1974 and August - October 1975. Numerals give number of data 
points at each position. Best-fit linear regression line is shown. Vegetative 
cover was estimated at 40 % . Dashed line encloses anomalous data 
points which were investigated more fully. 
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FRACTIONAL SATURATION OF TOP ZONE 


FIGURE 12. Daytime brightness temperature versus fractional saturation of top zone 
for five stations in the western part of the United States corn belt for 
September - October 1974 and August - October 1975. Numerals give 
number of data points at each position. Best-fit linear regression line is 
shown. Vegetative cover was estimated at about 55 % . 
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FIGURE 13. Daytime brightness temperature versus fractional saturation of top zone 
for three stations in the southern part of the United States winter wheat 
belt for September - October 1974 and August - October 1975. Numerals 
give number of data points at each position. Best-fit linear regression line 
is shown. Vegetative cover was estimated at about 40% . 
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FIGURE H: Andlyzecl ESMR 6 map of T^or over MinnebOia, 9 Dec 1975, 1718 GMT. Numbers in circles 
are snow depth reports in centimeters 
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